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Introduction 

The  inositol  signaling  enzyme,  phospholipase  C-yl,  (PLC-yl),  catalyses  the 
hydrolysis  of  phosphatidylinositol  4,5-bisphosphate  (PIP2 ),  to  the  two  second 
messengers,  inositol  3,4,5-trisphosphate  (IP3)  and  diacylglycerol  (DAG).  These  two 
second  messenger  play  important  roles  in  the  regulation  of  cellular  function.  Abnormal 
levels  of  DAG  have  been  associated  with  tumor  formation  and  growth.  Overexpressed 
and  activated  of  PLCyl  has  been  seen  in  breast  cancer  tumors  (1,2).  One  of  the  results 
of  overexpression  and  activation  of  PLCyl  in  breast  tumors  appears  to  be  increased  cell 
migration.  Therefore,  it  has  been  proposed  that  PLCyl  may  play  a  role  in  tumor 
metastasis  (3).  Inhibition  of  PLCyl  function  by  various  agents  has  been  shown  to  prevent 
growth  factor-induced  cell  migration  (3),  which  makes  this  enzyme  an  attractive  target 
for  novel  anti-cancer  therapies.  Since  PLCyl  also  plays  a  role  in  normal  cellular 
function,  the  global  inhibition  of  this  enzyme  is  not  desirable.  Also,  PLCyl  is  one 
member  of  a  family  of  isozymes  that  catalyze  the  hydrolysis  of  PIP2.  The  other  isozymes 
have  been  implicated  in  playing  roles  in  the  cell  different  from  PLCyl.  However,  all  of 
the  iosyzmes  in  this  family  are  believed  to  have  the  same  catalytic  mechanism,  only  the 
modes  of  regulation  are  different.  Any  catalytic  inhibitor  of  PLCyl  will  also  inhibit  the 
other  isozymes.  Therefore,  agents  designed  to  inhibit  the  regulation  of  PLCyl  would  be 
ideal.  The  purpose  of  this  Concept  grant  was  two-fold;  1)  to  investigate  if  kinetic 
simulation  programs  based  on  the  current  kinetic  model  for  the  hydrolysis  of  PIP2  by 
PLCyl  using  existing  experimental  data  can  simulate  the  kinetic  behavior  of  PLCyl  and 
2)  to  use  the  kinetic  simulation  program  to  identify  constriction  points  for  the  production 
ofPIP2. 

Body 

Statement  of  Problem 

In  response  to  EGF,  PLCyl  undergoes  translocation  from  the  cytosol  to  the 
membrane  and  becomes  activated.  Though  calcium  is  required  for  activity,  the 
translocation  process  does  not  appear  to  calcium  dependent.  Review  of  the  literature 
indicates  that  levels  of  certain  anionic  lipids  maybe  the  regulator  of  the  membrane 
binding  step.  In  addition  to  be  possible  regulators  of  membrane  binding,  anionic  lipids 
also  appear  to  regulate  the  production  of  IP3. 

It  is  not  clear  how  these  calcium  or  anionic  lipids  regulate  PLCyl  activity  nor  is  it 
clear  which  reaction  steps  are  being  regulated  by  these  agents.  The  first  step  in  being 
able  to  design  agents  to  regulate  the  amount  of  PIP2  hydrolyzed  by  PLCyl  is  to  be  able  to 
identify  the  steps  that  control  the  flux  through  the  pathway.  Using  current  experimental 
systems,  it  is  not  possible  to  determine  if  which  steps  in  PIP2  hydrolysis  by  PLCyl  are 
responsible  for  controlling  flux  through  the  system 

Background  and  Significance 

The  current  kinetic  model  for  PIP2  hydrolysis  by  PLCyl  was  proposed  based  on  in 
vitro  kinetic  experiments  (Fig  1)  (4).  In  this  model  cytosolically  located  PLCyl  (E)  first 
binds  to  the  membrane  (M)  to  form  the  enzyme  bound  complex  (EM)  (step  1).  Once  on 
the  membrane  surface,  substrate  (S)  can  diffuse  into  the  active  site  after  which  hydrolysis 
of  PIP2  occurs  (step  2).  Successive  cycles  of  diffusion  of  PIP2  into  the  active  site  with 
subsequent  PIP2  hydrolysis  can  occur  before  the  enzyme  dissociates  from  the  membrane 
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surface.  This  type  of  reaction  can  be  regulated  at  either  step  1,  membrane  binding  or  step 
2,  substrate  catalysis.  Experiments  conducted  so  far  has  demonstrated  little  regulation  of 
step  1,  which  is  consistent  with  the  literature. 


Figure  1.  Model  of  production  of  IP3  by  PLCyl.  Step  1;  Cytosolic  PLCyl  (Ec)  associates 
with  the  membrane  to  form  the  enzyme: membrane  complex  (EM).  Step  2;  The 
enzyme:membrane  complex  (EM)  hydrolysis  PIP2  (S)  to  IP3  (P). 

In  the  simplest  form,  this  reaction  would  follow  an  order,  Bi  Bi  kinetic 
mechanism.  Since  the  membrane  is  assume  to  be  saturating,  this  reaction  would  simplify 
to  Michaelis-Menton  reaction  mechanism,  generating  a  hyperbolic  velocity  curve  in 
response  to  increasing  substrate. 

1 

Eqn(l)  E  +  M  =  EM  +  S  - ►EM  +  P  V  0-5 


0 

s 

- 


However,  experimental  data  measuring  the  amount  of  IP3  formed  by  PLCyl  as  a  function 
of  PIP2  concentration  generated  sigmodal  curves  (4).  Indicating  that  catalysis  of  PIP2  by 
PLCyl  follows  a  mechanism  that  involves  binding  of  an  activator  (A)  to  a  site  separate 
from  the  active  site.  This  suggests  that  step  2,  the  hydrolysis  of  substrate,  is  the  regulated 
step.  Further  support  for  the  regulation  of  step  2  by  an  activator  was  obtained  when  it 
was  shown  that  other  anionic  lipids,  lipids  that  can  not  be  hydrolyzed  by  PLCyl,  could 
increase  PIP2  hydrolysis.  In  the  presence  of  the  anionic  lipids,  the  reaction  kinetics 
switched  from  a  sigmodal  function  to  hyperbolic  function.  The  parameter  affected 
appeared  to  be  substrate  binding,  in  the  presence  of  activator  the  half-maximal  binding  of 
PLCyl  decreased  10-fold.  The  change  in  the  kinetic  behavior  of  PLCyl  in  the  presence  of 
the  lipid  activator  suggests  that  the  enzyme  exists  in  at  least  two  conformations  on  the 
membrane  surface,  inactive  and  active.  As  a  result  of  these  experiments  equation  1  was 
written  to  include  an  inactive  and  active  form  of  the  PLCyl  on  the  membrane  surface. 

Eqn  (2)  E  +  M  =  EM  +  A  =  EMA  +  A  =  EMA2  +  A"  =  EMAn  +  S  =  EMA"  +  P 

However,  these  kinetic  experiments  did  not  directly  measure  the  membrane 
binding  step.  Experiments  designed  to  measure  the  membrane  binding  step,  found  that 
membrane  composition  did  not  affect  membrane  binding,  but  did  affect  calcium  binding. 


Data  showed  that  PLCyl  had  two  single  ion  binding  sites  for  calcium  (Miller  and  Jones, 
unpublished  results).  Site  1  had  an  low  affinity  for  calcium  (0.02  pM)  and  is  presumed  to 
be  the  active  site,  Site  2  was  found  to  have  an  variable  affinity  for  calcium.  In  solution, 
the  calcium  affinity  was  10  pM,  in  the  presence  of  neutral  phospholipids  the  affinity  was 
1  pM  and  in  the  presence  of  activating  lipids  the  affinity  was  0.04  pM  (Jones, 
unpublished  results).  These  data  gave  experimental  evidence  to  the  conclusions  from  the 
kinetic  experiments  that  there  are  at  least  two  conformations  of  membrane  PLCyl , 
depending  on  the  presence  of  activating  lipid.  However,  these  experiments  can  not  allow 
us  to  conclude  if  the  change  in  calcium  affinity  in  the  presence  of  activating  lipid  will 
effect  the  production  of  IP3, 

Method 

Metabolic  control  analysis  (MCA)  is  designed  to  aid  in  the  solving  of  these  type 
of  problems.  Using  user  defined  parameters,  MCA  will  allow  the  determination  of  the 
flux  through  a  system  at  steady-state.  Analysis  of  perturbations  at  various  steps  will 
allow  the  determination  of  which  steps  are  responsible  for  controlling  the  flux  through 
the  system.  Currently  an  number  of  freeware  programs  are  available.  In  the  proposal, 
simulations  using  the  current  kinetic  mechanism  and  experimental  derived  kinetic 
parameters  for  the  hydrolysis  of  PIP2  by  PLCyl  were  run  using  the  freeware  Gepasi. 

From  these  simulations  the  followed  questions  were  asked: 

•  Do  the  simulations  resemble  experimental  obtained  results? 

•  What  step  most  influence  the  amount  of  PIP2  hydrolyzed  ? 

Results 

For  the  initial  simulations,  each  reaction  step  in  the  hydrolysis  of  PIP2  by  PLCyl 
was  entered  in  MCA  program.  All  competing  reaction  steps  were  included,  This 
situation  likely  mimics  the  actual  events  (Fig.  2).  The  kinetic  parameters  used  were  those 
that  had  be  previously  determined  experimentally  (Table  1).  A  kinetic  mechanism  based 
on  the  experimental  results  was  selected  for  each  step.  The  results  of  these  simulations 
did  not  mimic  the  experimental  data.  Further  analysis  revealed  that  using  the  reaction 
scheme  as  written,  the  program  could  not  determine  a  steady-state  for  the  reaction 
scheme.  Breaking  the  reaction  scheme  into  sub-reactions  revealed  that  the  problem 
appeared  to  lie  with  formation  of  product  from  substrate.  Attempts  to  solve  this  problem, 
have  so  far  been  unsuccessful.  It  is  possible  that  kinetic  mechanisms  propose  for  each 
step  in  the  reaction  scheme  is  incorrect.  Selecting  different  mechanisms  did  not  solve  the 
problem.  It  is  also  possible  the  kinetic  parameters  derived  are  incorrect  The  MCA 
program  contains  functions  that  will  allow  the  determination  of  kinetic  parameters  from 
the  fitting  of  experimental  data.  Currently,  all  previously  collected  experimental  data  is 
being  converted  into  a  format  that  can  be  read  by  the  program.  Once  completed,  the  data 
will  be  fitted  using  the  MCA  program.  It  is  also  possible  that  more  reaction  steps  are 
required  for  the  hydrolysis  of  PIP2  than  have  been  proposed.  If  this  is  the  case  addition 
experiments  will  need  to  be  conducted. 
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Table  1.  Kinetic  mechanisms  and  parameters  used  for  reaction  scheme  depicted  in 


Reaction 

1 - ; - ; . . *  _ 

Kinetic  Mechanism 

Kinetic  parameters 

E  <->  Eca 

Mass  action 

Ks-  4  fiM  (Jones,  unpublished 
results) 

Eca  <->  EMCa 

Mass  action 

Ks  =  250  pM  (Jones,  unpublished 
results) 

EMCa<->  ESCa 

Catalytic  activator 

Km  =  0.2  mf,  Kc  =  1  mM  (4) 

E  <->  ExM 

Mass  action 

Ks  -  250  pM  (Jones,  unpublished 
results) 

ExM  <->  ExMCa 

Mass  action 

Kg  =  1  pM  (Jones,  unpublished 
results) 

ExMCa  <->  ExACa 

Catalytic  activator 

Km  =  0.2  mf,  Ks  =  100  pM  (4) 

ExM  <->  ExA 

Hyperbolic  modifier 

Ks  =  100  pM  (4) 

ExA  <->  ExAxCa 

Mass  action 

K,  =  0.4  uM  (Jones,  unpublished 
results) 

S  <->  P 

Explicit  Michaelis-  Menton 

Km  =  0.05  mf  (4) 

Fig  2.  Possible  reaction  scheme  describing  the  formation  of  IP3  by  PLCyl.  Starting  with 
PLCyl  located  in  the  cytosol  (E),  depicted  here  are  possible  reaction  schemes  for  the 
formation  of  1P3.  This  reaction  schemes  are  based  on  the  available  experimental  data. 
PLCyl  (E);  membrane  (M);  calcium  (Ca);  A  (anionic  lipid  activator);  PIP2  (S);  IP3  (P). 


It  was  possible  to  model  the  membrane  binding  steps  of  the  reaction  scheme 
(Figure  3).  The  simulations  run  mimicked  data  that  was  obtained  experimentally. 
Following  is  an  examples  of  the  type  of  questions  that  can  be  answered  from  these  tvpe 
of  simulation  studies.  In  the  cell,  activation  of  enzyme  cascades  by  calcium  occur  around 
0.3  [iM.  From  a  simulations  run  under  conditions  described  in  Table  2,  it  can  be  seen 
that  significant  concentrations  of  active  complex  are  only  formed  when  membrane 
binding  precedes  calcium  binding  to  site  2  (>80%  for  membrane  binding  first,  in  the 
absent  of  anionic  lipids,  as  compared  to  <  5%  for  calcium  binding  first  (Figure  4)  and 
>90%  for  membrane  binding  in  the  presence  of  activating  lipid,  as  compared  to  <  1%  for 
calcium  binding  first  (Figure  5)).  More  interesting  is  the  simulation  that  demonstrate 
where  half-maximal  formation  of  active  enzyme  complex  occurs  (Figure  6  and  7).  For 
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reaction  1,  half-maximal  formation  of  active  complex  occurred  at  3-4  pM  (EMCa)  and 
for  reaction  2,  the  half -maximal  formation  of  the  enzyme  complex  occurred  at  luM  in  the 
absence  of  activating  lipid  (Figure  6)  and  0.3  pM  in  the  presence  of  activating  lipid 
(Figure  7).  The  flux  due  to  the  formation  of  active  complex  was  compared  for  the  system 
(Figure  8  and  9).  In  the  presence  of  activating  lipid,  large  changes  in  flux  occurred  at 
0.075  pM  and  0.2  pM  (Figure  9),  whereas  in  the  absence  of  activating  lipid,  significant 
change  is  the  flux  occurs  at  1  pM  (Figure  7).  These  simulations  indicate  that  in  the 
presence  of  activating  lipid  the  formation  of  active  PLCyl  is  sensitive  to  changes  in 
calcium  concentration  and  furthermore,  that  changes  in  calcium  concentration  within  the 
physiological  range  could  have  a  profound  effect  on  active  complex  formation.  However 
the  problem  with  these  simulations  is  that  under  saturating  concentrations  of  membrane, 
such  as  would  be  found  in  the  cell,  most  of  the  enzyme  would  be  bound  to  the  membrane. 
This  is  not  the  case  in  cells.  Under  resting  conditions  <20%  of  PLCyl  is  associated  with 
the  membrane  (5).  The  discrepancy  between  what  occurs  in  the  cell  and  the  simulation 
results  obtained  from  in  vitro  experiments  may  indicate  that  steps  are  missing  in  the 
reaction  scheme.  One  possible  missing  step  is  that  PLCyl  prefers  to  bind  to  lipid 
domains  (6).  Kinetic  reactions  are  being  written  to  test  this  possibility. 

Table  2.  Kinetic  Mechanism  and  parameters  for  the  reaction  scheme  depicted  in 


Fig- 3. 


Reaction 

Kinetic  Mechanism 

Kinetic  Parameters 

E  <->  Eca 

Mass  action 

K,  -  4  pM  (Jones, 
unpublished  results) 

Eca  <->  EMCa 

Mass  action 

Ks  =  250  pM  (Jones, 
unpublished  results) 

E  <->  ExM 

Mass  action 

j 

K8  =  250  pM  (Jones, 
unpublished  results) 

ExM  <->  ExMCa 

Mass  action 

Ks=  0.4  pM  (Jones, 
unpublished  results) 

+Ca 


reaction  2 


reaction  1 


Figure  3.  Competing  reactions  for  the  formation  of  the  active  complex  of 
PLCyl.  The  cytosolic  located  enzyme  (E)  can  either  bind  calcium  (Ca)  in  the 
cytosol,  then  associate  with  the  membrane  (M)  to  form  the  active  complex 
(EMCa)  [reaction  1]  or  associate  with  the  membrane  then  bind  calcium  to  form 
the  active  complex  ExMCa  [reaction  2]. 
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These  simulations  suggest  focusing  experimental  efforts  on  the  binding  of  anionic 
lipid  to  PLCyl  would  be  beneficial.  The  simulations  indicate  that  by  inhibiting  the 
anionic  lipid  induce  changes  in  calcium  affinity,  would  regulate  the  flux  through  the 
system  by  shifting  the  formation  of  active  enzyme  complex  to  calcium  concentrations 
greater  than  that  physiologically  seen. 

In  addition  to  the  simulations  presented  in  this  report,  simulations  conducted 
under  the  following  conditions  are  ongoing: 

1 .  What  effect  increases  in  membrane  binding,  like  observed  with  tyrosine 
phosphorylated  PLCyl,  affects  membrane  binding  and  formation  of  active  complex 
enzyme. 

2.  What  effect  increases  in  tyrosine  phosphorylated  enzyme  in  response  to  EGF  has  on 
membrane  binding  and  formation  of  active  enzyme  complex. 

3.  What  effect  formation  of  lipid  clusters,  such  a  rafts  formed  in  cells,  would  have  on 
membrane  binding. 

To  date  reactions  have  been  written  for  all  three  and  simulations  run  for  1  and  2,  but  the 
reactions  have  not  been  linked  together  yet. 
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Key  Research  Accomplishments 

The  key  research  accomplishments  for  this  projected  were: 

•  Various  MCA  programs  were  evaluated  and  one  was  found  that  was 
suitable  for  the  PLCyl  analysis  (Gepasis). 

•  Reaction  equations  for  the  binding  the  PLCyl  binding  step  and 
enzyme  activation  were  written  and  evaluated. 

•  Simulations,  using  kinetic  equations  from  experimental  data,  of  the 
PLCyl  binding  and  activation  step  mimicked  the  experimental  results. 

•  Simulations  of  the  PLCyl  binding  and  enzyme  activation  indicate  that 
membrane  binding  probably  occurs  before  calcium  binding. 

•  Simulations  of  the  PLCyl  binding  and  enzyme  activation  indicate  that 
the  most  probable  control  point  for  this  reaction  is  the  calcium  binding 
step. 

•  Inhibition  of  the  calcium  binding  step  either  by  controlling  calcium 
concentration  or  decreasing  the  affinity  of  calcium  for  PLCyl,  did  not 
affect  membrane  binding,  but  did  effect  formation  of  the  active 
enzyme. 

Reportable  Outcomes 

To  date:  none 

Conclusions  and  Future  directions 

Simulations  conducted  so  far  has  shown  that  this  approach  could  be  quite  useful 
in  the  identification  of  possible  regulatory  steps  involved  in  the  EGF-induced  hydrolysis 
of  PIP2  by  PLCyl.  Based  on  the  results  of  the  membrane  binding  simulations,  additional 
experiments  will  be  conducted  to  investigate  the  change  in  calcium  binding  affinity  in  the 
presence  of  anionic  lipids.  Our  first  attempt  to  model  the  reaction  scheme  from 
membrane  binding  to  formation  of  product  was  not  successful.  The  best  approach  will  be 
to  break  down  the  reaction  scheme  into  subunit  that  correspond  to  the  experimental  data. 
Once  the  subunit  reactions  have  been  modeled  and  the  results  of  the  simulation  resemble 
the  experimental  data,  then  the  subunit  reactions  can  be  linked  together. 

The  advantages  to  being  able  to  simulate  EGF-induced  hydrolysis  of  PIP2,  is  that 
the  simulation  will  take  into  account  all  known  perturbations  on  the  system.  Experiments 
conducted  to  data  have  looked  at  only  certain  parts  of  the  reaction  scheme  under 
controlled  conditions.  Regulatory  steps  of  an  isolated  system  may  not  be  relevant  when 
all  perturbations  to  a  reaction  are  taken  into  account.  The  use  of  MCA  may  indicate  novel 
regulatory  steps  that  can  be  targeted  for  rational  drug  design.  The  disadvantage  of  MCA 
is  that  it  is  a  computer  simulation  and  may  have  no  relevancy  to  what  occurs  in  the  cell. 

It  will  be  important  to  test  the  results  of  the  simulation  in  both  isolated  controlled  systems 
and  in  cells. 


16 


* 


References 

1.  Arteaga,  C.L.  et.  al  (1991)  Elevated  content  of  the  tyrosine  kinase  substrate 
phospholipase  Cyl  in  primary  human  breast  carcinomas.  Proc.  Natl.  Acad  Set. 
88:10435. 

2.  Sonderquist,  A.M.,  (1992)  Elevataed  membrane  association  of  phospholipase  C-yl  in 
MDA  468  mammary  tumor  cells.  Cancer  Res.  52:4526 

3.  Kassis,  J.  Moellinger,  J.,  Lo,  H,  Greenberg,  N.M.,  Kim,  H.-G.,  and  Wells,  A.  (1999) 
A  role  for  phospholipase  Cy-mediated  signaling  in  tumor  cell  invasion.  Clin.  Ca.n 
Res.  5:2251. 

4.  Jones,  G.  A.,  and  Carpenter,  G.C.,  (1993)  The  Regulation  of  Phospholipase  C- 
gammal  by  Phosphatidic  Acid,  J.  Biol.  Chem.  268:20845. 

5.  Todderud,  G.,  Wahl,  M.,  Rhee,  S.  G.,  and  Carpenter,  G.  C.,  (1990)  Stimulation  of 
Phospholipase  C-yl  Membrane  Association  by  Epidermal  Growth  Factor,  Science 
249:296. 

6.  Pike,  L.  J.,  and  Miller,  J.  M.  (1998)  Cholesterol  depletion  delocalizes 
phosphatidylinositol  bisphosphate  and  inhibits  hormone-stimulated 
phosphatidylinositol  turnover,  J.  Biol.  Chem.  273:22298. 


17 


